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ABSTRACT 

We aim to understand the correlation between cloud formation and alkali line formation in 
substellar atmospheres. We perform line profile calculations for Na I and K 1 based on the 
coupling of our kinetic model for the formation and composition of dust grains with ID radia- 
tive transfer calculations in atmosphere models for brown dwarfs and giant gas planets. The 
Na 1 and K 1 line profiles sensibly depend on the way clouds are treated in substellar atmo- 
sphere simulations. The kinetic dust formation model results in the highest pseudo-continuum 
compared to the limiting cases. 
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1 INTRODUCTION 

I Dust influences the environment from which it forms by consum- 
, ing elements from the gas phase. Hence, dust selectively alters the 
' local metallicity in the atmosphere of a brown dwarf as well as in 
planetary atmospheres (Woitke & Helling 2004, Helling, Woitke 
■ & Thi 2008). Dust clouds furthermore efficiently absorb photons 
' which can result in /jm-broad absorption features (Helling et al. 
', 2006, 2007; Gushing et al. 2006, Burgasser et al. 2007). The ab- 
sorbed energy is redistributed inside the lattice structure of the solid 
grain causing an isotropic irradiation of IR photons into the atmo- 
, sphere, a process, which considerably alters the atmospheric tem- 
perature profile (Tsuji et al. 1996, Tsuji 2005, Ackerman & Marley 
2001, AUard et al. 2001). Both processes, the selective alteration 
of local element abundances and the changing atmospheric temper- 
ature, have a distinct influence on the gas-phase chemistry which 
determines the atomic and molecular concentrations. These intri- 
cate effects are important even for elements not or only marginally 
involved into the dust formation process as we will show in this 
letter. 

Alkali lines are dominant features in the optical spectrum of 
brown dwarfs and the far Une wings of the most abundant alkali 
metals, sodium and potassium, determine the pseudo-continuum in 
this spectral range (Johnas et al. 2007a, AUard et al. 2007). Another 
alkali, Li, is used to determine fundamental parameter of substellar 
objects (e.g. Martin, Rebolo & Maguzza 1994, Johnas et al. 2007c). 
They all can be a powerful tool to study the cloud formation in sub- 
stellar objects if combined with an atmosphere simulation which 
allows to correlate the line profiles with an analysis of the atmo- 



sphere regarding the Pgas, v^^m) structure {T - gas temperature, 
Pgas - gas pressure, Vconv - convective velocity) and the chemistry. 

In this letter we show that the modelling of dust formation 
has a strong impact on the line shapes and the depth of the Na I 
and K I alkali lines in dust-enshrouded L-dwarfs. We utilise our 
kinetic model for dust cloud formation in atmosphere simulations 
of substellar objects (brown dwarfs and extrasolar planets) in com- 
parison to two limiting cases in order to demonstrate the sensitive 
dependence of alkali line profiles on the local {T, pg^s, e,) profile (e, 
- element abundance) of the atmosphere. We show that the treat- 
ment of dust cloud formation has a much stronger impact on the 
alkali line profiles than the details of the line profile calculation. 



2 ATMOSPHERE MODELS 

We present our studies based on simulations carried out with the 
general purpose model atmosphere code PHOENI^Q, version 15, 
(Hauschildt & Baron 1999). We combine the study of line pro- 
files of Johnas (2007) with the kinetic dust cloud modelling by 
(Helling & Woitke 2007) which has recently been incorporated 
into PHOENIX by Dehn (2007) (also Dehn et al. 2008). We con- 
centrate here on the following two different line profile approaches 
for which the model atmospheres are calculated: 

• IMPACT - LINE profiles: 

van der Waals profiles applying the impact approximation for near 
line wings (Schweitzer et al. 1996) 
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This is the same code as used by Allard et al. (2001) except for the dust 
chemistry and line-profile calculations. 
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• modernI - LINE profiles: 

Detailed non-analytical, semi-classical line profiles (Allard et al. 
2005, Allard & Spiegelman 2006, Johnas et al. 2006) including H2 
in two symmetries, C2, and C„,., and He as perturber. The alkali line 
profiles are represented by two terms, one describing the near line 
wing with the impact approximation (Baranger 1958b, Royer 1971) 
and the other describing the far line wing with the one-perturber 
approximation in the density expansion (Royer 1971).|^ 

For our study of the influence of different dust model approaches 
on alkali line profiles, we compare simulations for the two limit- 
ing cases CoND and Dusty (Allard et al. 2001), and our new Drift 
module (Dehn et al. 2008) which includes a detailed modelling of 
dust formation: 

• CoND approach: 

The CoND approximation assumes the dust to be in chemical and 
phase equilibrium with the surrounding gas phase. The dust is ig- 
nored as opacity source because it is assumed that after the dust 
grains have formed they sink down below the photosphere, thus re- 
sulting in a depleted gas phase with no dust absorption or emission. 
This scenario was used to model T-type brown dwarfs (e.g. Lodieu 
et al. 2007) or planetary objects (e.g. Seifahrt et al. 2007) with Tejf 
lower than ~1400K but does not allow to consider the cloud forma- 
tion process. These models are not realistic for the parameter range 
considered here and are included for comparison only. 

• Dusty approach: 

The dust is assumed to have formed in chemical and phase equilib- 
rium, as in the Cond models but is assumed to remain at the place of 
formation. In both Cond and Dusty, the amount of elements bound 
in grains is subtracted from the gas phase. No dust settling is taken 
into account. The interstellar grain size distribution is used for the 
dust opacity calculation. Dusty models were used at higher Ted rep- 
resenting the case of a thick cloud layers inside the atmosphere of 
L- dwarfs (e.g. Leggett et al. 2001) and do not consider the actual 
formation process. 

• Drift approach: 

The Drift package (Dehn et al. 2008) is based on a kinetic treat- 
ment of dust formation (Woitke & Helling 2003, 2004; Helling & 
Woitke 2006; Helling, Woitke & Thi 2008). A stationary dust for- 
mation process is assumed for application in static model atmo- 
spheres. In these model atmospheres, seed particles (here: Ti02) 
nucleate from the gas phase if an appropriate super-saturation is 
achieved, subsequently grow a mantle made of various compounds 
(Si02[s], AI2O3M, Fe[s], MgO[s], MgSiOsls], Mg2Si04[s], and 
Ti02[s]) by 32 chemical surface reactions (Dehn 2007), settle grav- 
itationally, and evaporate. The dust formation cycle is completed by 
convective overshooting of uncondensed material which is mod- 
elled by an exponential decrease of the mass exchange frequency 
into the radiative zone (Ludwig et al. 2006). The material composi- 
tion of the grains, the amount and the size of dust particles formed 
- which are needed to evaluate the cloud's opacity in the radiative 
transfer calculation - are results of our model. 



3 RESULTS 

We demonstrate our results for a L-type brown dwarf with T^g = 
2000 K, log(g) = 5.0, and solar composition (Anders & Grevesse 

^ Johnas et al. (2007a,b) show that alternative alkali near line wing pro- 
files (Mullamphy et al. 2007) result in negligible effects on the atmospheric 
structure and synthetic spectrum. 
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Figure 1. Left axis (solid lines) (T, Pgas) structures for the Drift-Phoenix, 
Dusty-Phoenix, and Cond-Phoenix models with the modem 1 alkali line pro- 
files for T^fi = 2000K, log(g) = 5.0. Right axis (dashed lines): Flux con- 
tribution function at 5880 A of the six different model (Drift, Dusty, and 
Cond models with impact and MODERNi line profiles). 
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Figure 2. The partial pressure p [bar] of Na and K as function of gas pres- 
sure pgas [bar] as results of the model runs for the Drift, Dusty, and Cond 
dust model approaches. 



1986). Figure [T] shows that the atmospheric (r, p)-structures are 
very similar for pgas > 1 bar for Drift-Phoenix and Dusty-Phoenix 
models, but can differ by up to 300K (i.e. 30 %) at lower pres- 
sures. All models agree well in the innermost atmospheric layers. 
The Drift-Phoenix temperature gradient is steeper compared to the 
Dusty-Phoenix case. The Drift-Phoenix model structure falls basi- 
cally between the two limiting cases Dusty and Cond. The different 
line profiles cause negligible differences in the (T, p)~structures for 
all three cloud models (Johnas 2007). 

In addition to the (J, p)-structures in Fig.[T] the flux contribu- 
tion function ^£fr (see Fuhrmeister et al. 2006) of each of the simu- 
lations is shown at 5880 A which is a line wing wavelength in the 
Nal D2 resonance line (Fig. [5). The maximum of is correlated 
to the location in the atmosphere at which the line wing is formed 
at that specific wavelength (Magain 1986, Fuhrmeister et al. 2006). 
Note, however, that '^f should only be considered for strong lines 
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Figure 3. Top: Drift-Phoenix model atmospheres for Teff = 2000K, log g 
=5. computed with the impact and modernI line profiles. Bottom: Relative 
changes of the top panel fluxes A(A, B) = (A ~ B)/A. 

like Nal and KI. Figure[T]demonstrates that in general the NaID2 
line wings at 5880 A form at comparable local gas pressures in 
the model atmospheres despite the different dust cloud approaches 
considered here. We furthermore find that Na and K appear with 
comparable concentrations in the line forming regime in all models 
(Fig.|2]l. This leaves the local temperature as only candidate for po- 
tential differences in the Na I and K I line profile for different dust 
cloud treatment. 

Figure [3] shows that the low resolution Drift-Phoenix spec- 
trum for two different line profile setups will basically differ in the 
optical and near-IR region which is, hence, the spectral region pre- 
dominantly influenced by the line profile treatment. Although one 
might have expected more effects in larger spectral ranges, the rel- 
ative changes shown in the bottom-panel of Fig. [3] are relatively 
small. On the contrary, the differences in the low-resolution spectra 
between the different dust approaches (Drift, Dusty) for fixed line 
profile treatment are quite remarkable (Fig.O. These strong dif- 
ferences are no surprise because the atmospheric (T, p)-structures 
diflFer substantially (Fig.[TJ, and therefore the amount of molecu- 
lar opacity carriers does change. For example, Dusty-Phoenix pro- 
duces a TiO-concentration of two orders of magnitude less than 
Drift-Phoenix at p„^^ x lO"'* bar. The temperature-dependent gas- 
phase composition does also influence the dust formation by pro- 
viding the constituents from which the dust forms. We find similar 
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Figure 4. Top: Dusty-Phoenix and Drift-Phoenix model atmospheres for 
Teff = 2000K, log g =5 computed with the modeml - line profiles. Bottom: 
Relative change of the top panel fluxes A(A,B) = (A - B)/A. 



discrepancies when comparing spectra between Drift-Phoenix and 
Cond-Phoenix models (not shown). 



3.1 Alkali lines 

The classical Dusty-Phoendc simulation result in a lower pseudo- 
continuum for both, the NalD and KI doublets compared to the 
Drift-Phoenix results (Fig. |5j. The Dusty-Phoenix models have a 
higher dust opacity than the Drift-Phoenix models which results 
in strongly differing (r, p)-structures. Hence, also the line shapes 
for the NalD and KI doublets differ considerably. It appears that 
the local gas pressure and the concentration of Na I and K I are 
very similar in the line forming region in these models but large 
differences occur in the local temperatures. Comparing Fig.[6]visu- 
alises that the Cond-Phoenix line profiles match up with the Drift 
results better in the far wings since here the differences in the local 
temperature are somewhat smaller (Fig. [T]l. Generally, the Drift- 
Phoenix atmosphere simulations, which are based on a detailed ki- 
netic cloud model, yield a stronger pseudo-continuum and deeper 
line cores. The similarity of the alkali partial pressures especially 
between the Drift and Dusty approach point to the local tempera- 
ture as main cause for the large differences in the line-profiles. All 
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Figure 5. Synthetic spectra displaying line profiles of the alkali species with the highest concentration in cool atmospheres, the Na I D (left) and K I (right) 
doublets, computed with Drift-Phoenix and Dusty-Phoenix models with impact - and modernl - line profiles. 
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Figure 6. The NalD doublet as computed with the impact- and modern 1 -line profile setups in Drift-Phoenix, Dusty-Phoenix, and Cond-Phoenix simulations. 
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figures show that the observable line shapes depend more on the 
dust model than on the line profile model. 



4 CONCLUSION 

We have demonstrated that the dust treatment in an atmosphere 
simulation for substellar objects has a large influence on the re- 
sulting line-profile of the alkalies Na I and K I. We conclude that 
the main cause for these large differences is the strongly changing 
local temperature which is directly linked to the dust treatment by 
the resulting dust opacity. The different treatments of the alkali line 
profiles result in negligible changes compared to the effect of dust 
treatment. 

The pseudo-continuum of alkali line-profiles are strongest if 
based on a microphysical dust model taking into account the kinetic 
nature of dust formation (Drift-Phoenix). This findings may foster 
studies on Li I and Rb I lines since Johnas (2007) has shown that 
in particular the far K I line wings can mask the Rb I doublet, but 
both the Na I and the K I overlap the Li I doublet cores. 
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